Diffusion of multi-isotopic chemical species in molten silicates

Diffusion experiments in a simplified Na2O–CaO–SiO2 liquid system are used to develop a general formulation for the fractionation of Ca isotopes during liquid-phase diffusion. Although chemical diffusion is a well-studied process, the mathematical description of the effects of diffusion on the separate isotopes of a chemical element is surprisingly underdeveloped and uncertain. Kinetic theory predicts a mass dependence on isotopic mobility, but it is unknown how this translates into a mass dependence on effective binary diffusion coefficients, or more generally, the chemical diffusion coefficients that are housed in a multicomponent diffusion matrix. Our experiments are designed to measure Ca mobility, effective binary diffusion coefficients, the multicomponent diffusion matrix, and the effects of chemical diffusion on Ca isotopes in a liquid of single composition. We carried out two chemical diffusion experiments and one self-diffusion experiment, all at 1250 C and 0.7 GPa and using a bulk composition for which other information is available from the literature. The self-diffusion experiment is used to determine the mobility of Ca in the absence of diffusive fluxes of other liquid components. The chemical diffusion experiments are designed to determine the effect on Ca isotope fractionation of changing the counter-diffusing component from fast-diffusing Na2O to slow-diffusing SiO2. When Na2O is the main counter-diffusing species, CaO diffusion is fast and larger Ca isotopic effects are generated. When SiO2 is the main counter-diffusing species, CaO diffusion is slow and smaller Ca isotopic effects are observed. In both experiments, the liquid is initially isotopically homogeneous, and during the experiment Ca isotopes become fractionated by diffusion. The results are used as a test of a new general expression for the diffusion of isotopes in a multicomponent liquid system that accounts for both self diffusion and the effects of counter-diffusing species. Our results show that (1) diffusive isotopic fractionations depend on the direction of diffusion in composition space, (2) diffusive isotopic fractionations scale with effective binary diffusion coefficient, as previously noted by Watkins et al. (2011), (3) self-diffusion is not decoupled from chemical diffusion, (4) self diffusion can be faster than or slower than chemical diffusion and (5) off-diagonal terms in the chemical diffusion matrix have isotopic mass-dependence. The results imply that relatively large isotopic fractionations can be generated by multicomponent diffusion even in the absence of large concentration gradients of the diffusing element. The new formulations for isotope diffusion can be tested with further experimentation and provide an improved framework for interpreting mass-dependent isotopic variations in natural liquids. 2014 Elsevier Ltd. All rights reserved. http://dx.doi.org/10.1016/j.gca.2014.04.039 0016-7037/ 2014 Elsevier Ltd. All rights reserved. ⇑ Corresponding author. Tel.: +1 5102255043. E-mail address: [email protected] (J.M. Watkins). 314 J.M. Watkins et al. /Geochimica et Cosmochimica Acta 139 (2014) 313–326